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Abstract 

Any-order pion inclusive distribution for a chaotic source in high energy 
collisions are given which can be used in both theory and experiment to 
analyze any-order pion interferometry. Multi-pion correlations effects on two- 
pion and three-pion interferometry are discussed. 
PACS number(s): 13.85 Hd, 05.30 Jp, 12.40 Ee. 

Hanbury-Brown and Twiss were the first who applied the Bose-Einstein (BE) corre- 
lations to measure the size of distant stars. The method was first applied to particle physics 
by Goldhaber et al.(GGLP) in 1959. Since then the size of the interaction region has 
been measured by numerous experiments in high energy collisions using different types of 
particles. Two-pion BE correlation is widely used in high energy collisions to provide the in- 
formation of the space-time structure, degree of coherence and dynamics of the region where 
the pions were produced || . Ultrarelativistic hadronic collisions provide the environment for 
creating dozens of pions |§-|7|], therefore one must take into account the effects of multi-pion 
BE correlations in those processes. Because the bosonic nature of the pion should affect 
the single and i-pion spectra and distort the i{i > 2)-pion correlation function, thus, it is 
very interesting to analyze the multi-pion BE correlation effects on i-pion interferometry 
P |15| . On the other hand, one also want to study higher-order pion interferometry directly 



to see what additional information can be extracted from higher-order pion interferometry. 
Now this aspect has aroused great interests among physicists |]o1,n jl^Ji6| HT8l . Unfortunately 



all present analyses of multi-pion correlations are based on pure "multi-pion interferometry 
formula" without considering the higher-order pion correlation effects on the lower-order 



pion interferometry |6U^,p^,p^-[T8| . Thus we urgently need new multi-pion interferometry 



formulas which can be used in both theory and experiment to analyze any-order pion in- 
terferometry. This is the main aim of this letter. In the letter, considering multi-pion BE 
correlations we derive new multi-pion correlations formulas which can be used to analyze 
any-order pion interferometry. Those new multi-pion correlations formulas are structurally 
similar to the previous formulas but with a modified source functions (See Eq. (21) for 
details). This warrants the validity of the formulas used in earlier studies of higher-order 
pion interferometry. Although we only study multi-pion BE correlations in this Letter, the 
results presented here are also held for kaon if the final state interactions are neglected. 
The general definition of the "pure" n pion correlation function C n (pi, • • •, p n ) is 

rt t n \ -fn(Pl) ■ " ") Pn) 

C n (Pl,- ■ -,Pn) = nB „, x , (1) 
lli=l -fllPiJ 

where P n (pi, • • -,p n ) is the probability of observing n pions with momenta {pi} all in the 
same n pion event. The n-pion momentum probability distribution P n (pi, • ■ -,p n ) can be 
expressed as PJfl|l 



Pn{Pl,' ■ -)Pn) = z2Pl,a(l)P2,a(2)---Pn,a(n), (2) 

(7 

with 

Pi tj = p(Pi, Pj) = / d 4 xg w (x, &*±M) e *(iH-wK (3) 



Here a(i) denotes the zth element of a permutation of the sequence 1, 2, 3, • • •, n, and the sum 
over a denotes the sum over all n! permutations of this sequence. g w (Y, k) can be explained 
as the probability of finding a pion at point Y with momentum k which is defined as JTTJ 



g w {Y, k)= J d A yf(Y + y/2)j{Y - y/2) exp(-zfc ■ y), (4) 

with 

J g w (x, k)d A xdk = n . (5) 
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Where j(x) is the current of the pion, no is the mean pion multiplicity without BE correlation 
mPA - From Eq. (1) and Eq. (2), the pure n-pion correlation functions can be expressed 



as 



re 



c n ( Pl ,-.-, Pn )=En— ■ w 

a j = l Pj,j 

The above correlation functions are widely used in both experiment and theory to ana- 
lyze multi-pion interferometry. But in high energy experiment, the pion multiplicity is so 
large that we must take into account the multi-pion correlations effects on lower-order pion 
interferometry. 

For nir events, considering the n pion correlations effect, the i-pion correlation function 



can be defined as 111 



where i^ n (pi, ■ ■ -,Pi) is the the normalized modified i-pion inclusive distribution in n pion 
events which can be expressed as 

P * (P1 ' ' ' ■' Pi) " /n^^p„( Pl ,---,p n ) ' (8) 



Now we define the function Gi(p, q) as [|T0|JTT 



Gi(p, q) = J p(p, pi)rfpip(pi, p 2 )dp 2 ■ • ■ p(Pi-2, Pi-i)rfpi-ip(Pi-i, q). (9) 

From the expression of P n (pi, • • -,Pn) (Eq. (2)), the one-pion to three-pion inclusive 
distribution can be expressed as fll 

1 1 n 

Pf(p) ttE%p)- w(n - •), (10) 



^ 1 n ' ^ 
P 2 n (Pi ; P2) = -7 TT^tHE C m (pi,pi) ■ Gi- m (p 2 ,p 2 ) 

+G m (pi, p 2 ) ■ G i _ m (p 2 , Pi)]w(n - i), (11) 
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J 1 n * — ^ * — m ~l 

^(Pi.fo.Ps) = -7- _ 7w- _ 9 , , s EtE E G Tn (pi,pi)-G fc (p a ,p a )-G < _ m ^ fc (ps,p 3 ) 
ft^n ij^i u^/ij i=3 m=1 fc=1 

+C m (pi, p 2 ) • G k (p 2 , Pi) ■ Gj_ m _ fc (p 3 , p 3 ) + G m (p 2 , p 3 ) • G k (p 3 , p 2 ) 
■Gi_ m _ fe (pi, pi) + G m (p 3 , Pi) ■ G fc (pi, p 3 ) ■ Gi_ m _ fc (P2, P2) + (12) 
G m (pi, p 2 ) ■ G fc (p 2 , p 3 ) • Gi^ m - k (p 3 , pi) + G m (pi, P3) 
■G k (ps, P2) ■ Gi_ m _ fe (p 2 , Pi)p(n - i), 

with 

1 r n 1 n r 

W M = — / II d Pk p n{Pi, ■ ■ ; Pn) = - E G(i)uj(n - i), C(i) = / Gj(p, p)dp. (13) 
n. j fc=1 n i=1 J 

Here a>(n) is the pion multiplicity distribution probability. 

Similar expression can be given for i(i < n) pion inclusive distribution. From the above 
method the i-pion correlation function can be calculated for n pion events. Experimentally, 
one usually mixes all events to analyze the two-pion and higher-order pion interferometry. 
Then the % pion correlation function can be expressed as [IS, 20] 

Ci (pij " " "i Pi) = n ^ iJVi(pj) . (14) 

Here the i-pion inclusive distribution, N{(pi, ■ ■ •, pi), can be expressed as 

N( , ggU u>(n) ■ n(n - 1) ■ ■ ■ (n - 1 + ( P i, • • pQ 

^(Pi> • • -,Pi) = ^ r~\ ' ( 15 ) 



with 



J N i (p 1 ,---,p i )f[dp i =<n(n-l)---(n-t + l)> . (16) 



Then the one-pion to three-pion inclusive distribution read: 



iV 1 ( P ) = E^(p, P ), (17) 



^2(Pi,Pa) = E G i(Pi'Pi)E G i(P2'P2) +E G i(Pi'P2)E G i(P2'Pi)' ( 18 ) 

i=l j=l i=l j=l 
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N 3 (pi, P2, Ps) = £ G i(Pi> Pi) £ G i(P2' Pa) £ G k(Ps, Ps) + 51 G i(Pi> P2) £ g j(P2> Pi) 

i=l j=l k=l i=l j=l 



00 00 00 00 



£ G k(P3, Ps) + £ G i(P3i Pi) I] G i(Pl> Ps) £ G k(P2, P2) + £ G;(p 2 , p 3 ) 

fc=l i=l j=l k=l i=l 

00 00 00 00 

I] Gj-(p 3 , P2) £ G k(Pl, Pi) + I] G i(Pl, P2) J! G i(P2> Ps) 

j=l k=l i=l j=l 

00 00 00 00 

£ G fc (p 3 , Pl ) + £ Cdh, p 3 ) £ G,(p3, Pa ) £ G fc (p 2 , Pi). (19) 
fe=l i=i j=i fc=i 

Similar expression for i(i > 3) pion inclusive distribution can be given. Following Ref. 

Il] , |Tol , PT|] , we define the following function 



H ij = H(p h p i ) = J2G k (p hPi ). (20) 

k=l 

Then the n pion inclusive distribution can be expressed as 

N n (Pl, • • -,Pn) = £ #1,7(1) #2,(7(2) * " * #n,<r(n)- (21) 
a 

Here a(i) denotes the ith element of a permutation of the sequence 1, 2,3, • • -,n, and the 
sum over a denotes the sum over all n\ permutations of this sequence. One of the interesting 
things about Eq. (21) is that it is very similar to Eq. (2). The only difference is that Eq. 
(2) only contains the the first term of ( p(pi, Pj) = Gi(pi,pj)). With the help of Eq. 
(21), the general form of z-pion correlation function (Eq. (14)) can be re-expressed as |16| 



Cf(pi,---,P,)=Ef[%^. ( 22 ) 
<7 j=l n i,o 

In the derivation of Eq.(22), we mention nothing about the structure of the source, so our 
results are in principle independent on the detail form of the source, i. e., one can use any 
kind source function (which may contain resonance and flow ) to study higher-order pion 
interferometry. Assumed that = \H^\exp{i(j)ij) , it is clear that two-pion interferometry 
does not depend on phase faj while which exists in higher order pion interferometry. So 



higher-order interferometry can be used to extract the information of phase [IT]. If the 
source distribution function is symmetric in the coordinate space we have 0^ = 0. Then 
multi-pion correlations contain the same information as two-pion interferometry does. In 



the following, we will use a simple model to study the multi-pion correlations effects on two- 
pion and three-pion interferometry. Similar to ref. P-Pf], we assume the source distribution 
function g(x,p) as 

9{x,p) = n ■ {^f' 2 exp{-r 2 /R 2 ){-^f' 2 exp{-p 2 /A 2 ). (23) 

Here R and A are parameters which represents the radius of the chaotic source and the 
momentum range of pions respectively. Using Eq.(9), Eq.(14) and Eq.(15), one can calculate 
any-order pion correlation function for any kind source distribution |22f . Csorgo and Zimayi 



have found an analytically solution for above special source distribution [|14i|15| which are 
quoted here: 

G„(pi,p 2 ) = h exp{-|[(7: /! pi " 1" /2 P2) 2 + (l! /2 P2 - l" /2 pi) 2 ]} (24) 
jn = nl^fl\K = ^^_ (25) 

with 

7± = ~{x± l) 2 ,x = RA. (26) 
Using this solution, one can calculat any-order pion interferometry for above source distri- 



bution. The analytical solution has the advantage over the previous method jlOlJlTl that it 
can be used in theory analyses. 

The two-pion interferometry for different mean multiplicity < n > is shown in Fig.l. 
It is clear that multi-pion correlation make the radius to become smaller. For three-pion 
interferometry, we choose variable Q 2 = (pi — P2) 2 + (P2 — P3) 2 + (P3 — Pi) 2 = Qi2 2 + q23 2 + 
q 31 2 and integrate the other eight variables, then we have three-pion correlation function 
C$(Q) as 

C <t>( \ S N 3(Pl, P2, Ps)S(Q 2 - qi2 2 ~ q23 2 ~ q31 2 )^Pl^P2^P3 

3lW /iV 1 (p 1 )iV 1 (p2)iV 1 (p3)5(Q 2 -qi2 2 -q23 2 -q3i 2 )rfpirfp2rfp3' 1 ' 

Three-pion interferometry for different (n) is shown in Fig.2. It is clear that as (n) becomes 

larger, the deviations of three-pion correlation from "pure three-pion "interferometry become 

larger. 



In the following, we will discuss the effects of resonance, flow and energy conservation on 
the above calculations. The effects of resonance on pure two-pion interferometry is not a new 
subject which has been extensive studied by different authors |[23|| . Effects of the resonance 
on higher-order pion interferometry based on pure multi-pion interferometry formulas are 
recently studied by Csorgo in Ref. Because Eq.(21) is similar to Eq.(2), all analyses 

in Ref. |16| can be applied here. In principle, resonance should not affected the multi-pion 



correlation formulas presented here, but due to the limited momentum resolution of the 
data we will have a modified n-pion BE correlation functions as presented in Ref. [|16| . The 



basic idea is that due to the limited resolution of the data, the contributions from the long 
lived resonance to the correlator are concentrated at lower relative momentum region which 
is not resoloved by two-pion interferometry. Thus the interference term (Gi(pi, pj), i ^ j) 
now mainly contains the contributions from directly emitted pions and pions from short 
lived resonances. While the single particle spectrum (Gi(pi,pi)) is not affected by the 
two-particle momentum resolution, so the intercept of n-pion interferometry is smaller than 
the ideal value n\. Based on pure two-pion interferometry formula, one has found that 
flow make the apparent radius derived from two-pion interferometry to become smaller 
23. It is expected that due to the multi-pion BE correlation effects, the apparent radius 



derived from two-pion interferometry will become more smaller. Considering the effects of 
the energy constraint on the multi-pion interferometry, I have derived a new multi-pion 
inclusive distribution according to the method presented in Ref. p5| . This new formula is 



similar to Eq.(21). If we integrated the energy from zero to infinity we will have Eq. (21) 
again. According to the definition of Eq.(14), in the calculation of Eq.(22), we actually 
mix all events (with different multiplicities and different energies). It means that in the 
experiment, one has already integrated all the possible energy of the pions, so it is not 
necessary for us to consider energy constraint effects here. 

Conclusion: In this letter, we have derived the i pion inclusive distribution and i pion 
correlation function for a chaotic source which can be used in experiment and theory to 
analyze multi-pion interferometry. For a simple model, multi-pion correlation effects on two- 



pion and three-pion interferometry were discussed. It was shown that for larger mean pion 
multiplicity, the deviations of three-pion and two-pion correlation from the "pure" three-pion 
and two-pion correlation became larger. The effects of resonance, flow and energy constraint 
on the multi-pion correlation formulas were discussed. 
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Figure Captions 

1. Multi-pion correlations effects on two-pion interferometry. The solid line corresponds 
to the result of pure two-pion interferometry. The dashed line and dotted line cor- 
respond to < n >= 20, 146 respectively. The input value of R and A are 5/m and 
0.25GeV. 

2. Multi-pion correlations effects on three-pion interferometry. The solid line corresponds 
to pure three-pion interferometry. The dashed line and dotted line correspond to 
< n >= 11,45 respectively. The input value of R and A are 3/m and 0.36GeV. 
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